<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"><html xmlns="http://www.w3.org/1999/xhtml"><head><title>R: An example pedigree data</title>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<link rel="stylesheet" type="text/css" href="R.css" />
</head><body>

<table width="100%" summary="page for l51"><tr><td>l51</td><td style="text-align: right;">R Documentation</td></tr></table>

<h2>An example pedigree data</h2>

<h3>Description</h3>

<p>The data contains data on 51 individuals in a pedigree. Below it is used for comparing results from various packages.
</p>


<h3>Usage</h3>

<pre>data(l51)</pre>


<h3>Format</h3>

<p>A data frame</p>


<h3>Source</h3>

<p>Morgan v3.</p>


<h3>References</h3>

<p>Morgan v3. <a href="http://www.stat.washington.edu/thompson/Genepi/MORGAN/Morgan.shtml">http://www.stat.washington.edu/thompson/Genepi/MORGAN/Morgan.shtml</a>
</p>


<h3>Examples</h3>

<pre>## Not run: 
km &lt;- kin.morgan(l51)
k2 &lt;- km$kin.matrix*2

# quantitative trait
library(regress)
r &lt;- regress(qt ~ 1, ~k2, data=l51)
names(r)
r
# qualitative trait
N &lt;- dim(l51)[1]
w &lt;- with(l51,quantile(qt,probs=0.75,na.rm=TRUE))
ped51 &lt;- within(l51, bt &lt;- ifelse(qt&lt;=w,0,1))
d &lt;- regress(bt ~ 1, ~k2, data=ped51)
d
# for other tests not shown here
set.seed(12345)
ped51 &lt;- within(ped51,{r &lt;- rnorm(N); bt[is.na(bt)] &lt;- 0})
library(foreign)
write.dta(ped51,"ped51.dta")

## End(Not run)
</pre>


</body></html>
